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c I ■ Abstract 



A model for the parton distributions of hadrons in impact parameter space has 
been constructed using soft gluon summation. This model incorporates the salient 
features of distributions obtained from the intrinsic transverse momentum behaviour 
of hadrons. Under the assumption that the intrinsic behaviour is dominated by soft 
gluon emission stimulated by the scattering process, the b-spectrum becomes softer 
and softer as the scattering energy increases. In minijet models for the inclusive 
cross-sections, this will counter the increase from a jet . 

The impact of parton scattering on the rise with energy of inclusive cross-sections was 
suggested by Cline and Halzen after such rise was first observed in proton-proton 
collisions at ISR . Minijet || |[ £| and eikonal minijet models H H, [7|, [U were subsequently 
developed to include an increasing number of partonic collisions in QCD resulting from the 
rapid rise in gluon densities. Recent measurements of photo- and hadro-production total 
cross-sections P, ITDf in energy regions where QCD processes dominate, confirmed the trend 



and were confronted with theoretical predictions obtaining varying degrees of success |TT . 
Typically, the eikonal (unitary) models for total hadron-hadron cross-sections are written 
as 

a total = 2j d 2 b[l - e -< b ^ 2 } (1) 
where the average number of collisions at impact parameter b is given by 

n{b,s) = A{b)a{s) (2) 
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In mini-jet models, a(s) = cr non - per t(s)+aj e t(s), with aj et to be calculated from perturbative 
QCD. A key ingredient of all models with a QCD component is the overlap function A(b), 
which describes matter distribution in impact parameter space. In some models ||, it 
has been assumed to be the Fourier transform of the product of hadronic form factors 



of the colliding particles. In other models [12], a gaussian shape has been assumed, thus 



relating A(b) to the intrinsic transverse momentum distribution of partons in the colliding 
hadrons. In either case, detailed information on A(b) relies on parameters to be determined 
case by case. Moreover, while direct measurements of the EM form factors are available for 
nucleons and pseudoscalar mesons, experimental information regarding photons or other 
hadrons such as vector mesons is lacking. The same observation applies to the intrinsic-pt 
interpretation for the spatial distribution of partons in vector mesons or photons. There 
is clearly need for a theoretical description of these distributions which would allow for 
a calculation of the mini-jet total cross-section (in the eikonal approximation) for photo- 
production and gamma-gamma collisions at future colliders. 

The aim of this paper is to provide a model for A(b) which could be applied to various 
cases of interest. We shall use Bloch-Nordsieck techniques to sum soft gluon transverse 
momentum distributions to all orders and compare our results with both the intrinsic 
transverse momentum approach as well as the form factor approach. In what follows, we 
shall first illustrate the Bloch-Nordsieck result and show that it gives a gaussian fall-off 
with an intrinsic transverse size consistent with Montecarlo models |12[] . We then calculate 



the relevant distributions and discuss their phenomenological application. 

We propose that in hadron-hadron collisions, the b-distribution of partons in the col- 
liding hadrons is the Fourier transform of the transverse momentum distribution resulting 
from soft gluon radiation emitted by quarks as the hadron breaks up because of the colli- 
sion. This distribution is obtained by summing soft gluons to all orders, with a technique 



amply discussed in the literature [O, The resulting expression is WB, p 



Fbm(K ± ) = i- / bdbJ (bK ± )e- h ^ M ^ (3) 



with 

2c F r M dki ,k 2 ^. M + x /M 2 -k 2 ± 



h(b; M, A) = / -±a s (-±)ln v t [l - J (k ± b)} (4) 

7T Jo k± A 2 M-JAP-k 2 , 



The choice of the hadronic scale M, which accounts for the maximum energy allowed in 
a single (k 2 = 0) gluon emission, plays a crucial role, just as it did for the Drell-Yan 
process, where the above expression has been successfully [17|, [H, HH, ^TJ used to describe 
the transverse momentum distribution of the time-like virtual photon or W-boson. In these 
cases[[l7|, |i~9| , the scale M was found to be energy dependent and to vary between \/Q^/4 
and y/Q 2 ~/2. We shall return to this determination later. 

The definition given in eqs.([L],|2]), requires for its consistency a normalized b-distribution, 

i.e. 

(5) 



J d 2 bA{b) 
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Our proposed Bloch-Nordsieck expression for the overlap function A(b), satisfying the 
above normalization, reads 

e -h(b;M,A) 

Abn = 27r/M6e-MWA) ( 6 ) 

An inspection of eq.([|), immediately poses the problem of extending the known asymptotic 
freedom expression for a s to the very small k± region. To avoid the small k±_ divergence in 
eq.(f|), it has been customary to introduce a lower cut-off in k± and freeze a s at k± = 0, 
i.e. to put 

(h 2 \ l2n 1 / 7 n 

as{±> 33-2N f \n[(k 2 ± + a 2 A 2 )/A 2 )] {> 



with a = 2 in ref. pi| . For applications where the scale M is large (e.g., W-transverse 
momentum distribution calculations) eq.(|4]) is dominated by the (asymptotic) logarithmic 
behaviour and the small fc^-limit, albeit theoretically crucial, is not very relevant phe- 
nomeno logically. However, this is not case in the present context, where we are dealing 
with soft gluon emission in \ow-p t physics (responsible for large cross-sections). The typ- 
ical scale of such peripheral interactions, is that of the hadronic masses, i.e. we expect 
M ~ 0(1 -T- 2 GeV) and the small k± limit plays a basic role. This can be appreciated on 
a qualitative basis, by considering the limit bM << 1 of eq. (^). In this region, we can 
approximate 1 — Jo(kb) pa b 2 k 2 /4, to obtain 

h(b; M, A) pa b 2 A (8) 

with 

We obtain a function h(b; M, A) with a gaussian fall-off as in models where A(b) is the 
Fourier transform of an intrinsic transverse momentum distribution of partons, i.e. 
exp(— kj_/4A 2 ). Note that the relevance of an integral similar to the one in eq.(^|) has been 
recently discussed in connection to hadronic event shapes |2"0 |. 



Our choice for the infrared behaviour of a s for a quantitative description of the dis- 
tribution in eq. does not follow eq. (|?]), but is inspired by the Richardson potential 
for quarkonium bound states p^| . In a number of related applications |2^, 23"| |, we have 



proposed to calculate the above integral using the following expression for a s : 

= (33 -I*,) ln[l +£(*(.)*] (10) 

which coincides with the usual one- loop expression for large (relative to A) values of k±, 
while going to a singular limit for small kj_. For the special case p = 1 such an a s coincides 
with one used in the Richardson potential |[24| , and which incorporates - in a compact 



expression - the high-momentum limit demanded by asymptotic freedom as well as linear 
quark confinement in the static limit. We have generalized Richardson's ansatz to values 
of p < 1. For 1/2 < p < 1, this corresponds to a confining potential rising less than linearly 
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with the interquark distance r. The range p ^ 1 has an important advantage, i.e., it allows 
the integration in eq.((|) to converge for all values of k±. A qualitative argument to justify 
the use of such (less singular) values for the parameter p follows. 

Assume a confining potential (in momentum space) given by "one-loop gluon" exchange 
term 

V(Q)=K(^f±), (11) 

where K is a constant calculable from the asymptotic form of a s . Let us choose the simple 
form 

" s(g2) = b ln[l + (QVA 2 ) p ]' (12) 
so that V(Q) for small Q goes as 

V(Q)^Q- 2 ^ 1+p \ (13) 
Eq.(|l3"D implies, for the potential, in coordinate space 



v{r) = I Wf elQrV ^ (14) 

as 

V( r ) -> (1/r) 3 • r (2+2p) ~ C r (2p - x) , (15) 

for large r (C is another constant). A simple check is that for p equal to zero, the usual 
Coulomb potential is regained. Notice that for a potenial rising with r, one needs p > 
1/2. A determination of p can be made following an argument due to Polyakov P7| . We 
minimize the potential, consistent with the constraint that the potential energy (or mass, 
M) increases as J^ 1//2 - ) for large J. ( A linear Regge trajectory means that the angular 
momentum J ~ M 2 ). The "confinement" energy plus the rotational energy ( all for large 
r and large J) reads 

V(r,J) ~{ J{J + l) ) + Cr^\ (16) 

Since the first term goes down with r and the second grows with r (for confinement, 
p > (1/2) is necessary), a minimum is guaranteed (for a fixed J). So we minimize the 
above potential, and find the position r = r a (J) where the minimum is, as a function of 
J, 

r Q (J) ~ J 2 /( 2 p +1 ), for large J. (17) 

Inserting this value of r Q , we find the minimum effective potential (or rest mass) as a 
function of J 

M^o)~J ( ^'. (18) 

A mass spectrum with linearly rising Regge trajectories, for which M 2 ~ J, implies the 
right side to grow like J^ 1//2 - ) , which gives 

p = (5/6). (19) 
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For the motivations given in [22] and repeated above, the value p = 5/6 was chosen 



in previous calculations of the transverse momentum distribution of Drell Yan pairs |22 



19| . Compared to eq.(^), we consider our present choice (i.e., eq.flTOD) physically more 
transparent as it leads to a direct quantitative estimate of soft-gluon generated intrinsic 
transverse momentum of partons in hadrons. With our ansatz, the intrinsic transverse 



Figure 1: Intrinsic transverse momentum of partons in a hadron for different hadron scales 
M, in units of A, for p=5/6. 

momentum becomes a calculable quantity rather than being an assumed one (see related 
discussion in Nakamura et al. [22 1). It can be obtained from the decomposition 

h(b; M, A) = hintrinsic + A(6; M, A) (20) 

with 



32 r A dk p . M + y/M 2 - k 2 



intrinsic—— TyrTT I T, /-, , m TT TTTg , o ~ ^0 ipk) ) , (21) 



33 - 2N f Jo k ln(l + p(g)J») M - VM 2 - k 2 
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where the contribution to the integral from the region k± < A determines the form and 
type of the "intrinsic transverse momentum" behaviour, as we discussed above. In this 
region, we now have 

b 2 2 

^intrinsic ~ ^ ^int (^2) 

with 

,9 32 1/1 , 2M\ a9 

< k > int = r — r + n — A 2 , (23) 

t mt 33 _2jv / (l-p) \2(l-p) A/ 

By comparison with eq.@, we see that A =< k 2 >i nt /4, which corresponds to an intrinsic 
transverse momentum of a few hundred MeV for A in the 100 MeV range and M < 1 GeV. 



For Nf = 3, we show in Fig.l the variation of y< k 2 >i nt as a function of M for a range 
of values of A. 

The "intrinsic" behaviour just discussed appears in the very small K± region, i.e. K± < 
A. In general, the full expression eq.(pO[) should be used. 



Having set up our formalism, we shall now examine its implications. The distribution 
A(b) depends upon the hadronic scale M in the function h(b). This scale depends upon 
the energy of the specific subprocess and, through this, upon the hadron scattering energy. 
In the calculation of the transverse momentum distribution of a lepton pair produced in 
quark- ant iquark annihilation |17[, the function e~ h ^ was the Fourier-transform of such 



distribution and the scale M was obtained as the maximum transverse momentum allowed 
by kinematics to a single gluon emitted by the initial qq pair of cm. energy \/§ in the 
process 

qq^g + i{Q 2 ) (24) 

Following ref. |I7|, for a lepton pair of mass squared Q 2 and rapidity y, this quantity is 
given by 



Js Q 2 1 
q ma x(s,y) = ^— ) = (25) 

2 s v 1 + zsinn z y 

with z = Q 2 /s. This would be the upper limit of integration in the function h(b) for each 
particular subenergy s of the quark-antiquark pair. The transverse momentum distribution 
can then be factorized from the overall hadronic cross-section, by substituting M in h(b) 
with the average of q max over all quark energies and configurations, with a weight 1/s, 
proportional to the cross-section. For a fixed value of the Feynman variable xp, one 
obtains 

Ul , I^I^S(x 1 -x 2 -x F )q maiB f(x 1 )f(x 2 ) 
M ^ XF) - S ^^5{x x -x 2 -x F )f{x 1 )f{x 2 ) W 

for given quark densities f(x). To understand the energy dependence of this scale, one can 
use a simple toy model, in which one approximates the valence quark densities with 1 / y/x 
and obtain from, eq.(|26|) at x F = 0, 

M(s, 0) « ^ ( 1 - 3^ Q«V~s (27) 
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which shows that M increases with energy towards the asymptotic value 2Q/3. This 
tendency was confirmed by numerical estimates with realistic quark densities J23|. An 
asymptotic increase of M with energy is also obtained for the case of interest here, i.e. 
parton-parton scattering contribution to the total cross-section. In the Drell-Yan case, one 
needed h(b) to calculate the transverse momentum distribution of the lepton pair, here 
we use it to evaluate the average number of partons in the overlap region of two colliding 
hadrons. In this case e~ h ^ is the F-transform of the transverse momentum distribution 
induced by initial state radiation in the process 

qq^ jet jet + X (28) 

where X can also include the quark- ant iquark pair which continues undetected after emis- 
sion of a gluon pair which stimulated the initial state bremsstrahlung. The jet pair in 



process (|28|) is the one produced through gluon-gluon or other parton-parton scattering 
with total jet-cross-section <jj et . We work in a no-recoil approximation, where the trans- 
verse momentum of the jet pair is balanced by the emitted soft gluons. Then the maximum 
transverse momentum allowed to a single gluon is still given by an expression similar to 
eq.(||), i.e. 

W(5) = ^(l-^) (29) 

except that now \fQ^ has been replaced by the jet-jet invariant mass ^Sj et} over which 
one needs to perform further integrations. In other words, this kinematic limit is the same 
as the one obtained for the gluon accompanying the Drell-Yan process mentioned above 
: there the quark-antiquark pair annihilates into a dilepton pair and a soft gluon, here 
it may continue on its way- undetected- after semi-hard emission of a parton pair which 
stimulated the soft gluon initial state bremsstrahlung. An improved eq.(|2|) now reads 

n(b, s) = n soft (b,s) + J2 J -^r J -^rfi( x i)fj( x 2) j dz J dp 2 t A BN {b } q max ) (30) 



where fi are the quark densities in the colliding hadrons, z = Sj et /(sxiX2), and is the 

differential cross-section for process ( |28l) for a given p t of the produced jets. In absence 
of a precise prescription of how to deal with the non-perturbative contribution to n(b,s), 
it is customary in these eikonal models to introduce a lower cut-off in the jet transverse 
momentum, usually indicated as ptmin to separate the mini-jet contribution from the soft 
part of the cross-section. 

Unlike the usual expressions for n(b, s), eq. ([!(]) does not exhibit factorization between 
the longitudinal and transverse degrees of freedom since the distribution Abn depends 
upon the quark subenergies. Factorization can be obtained however, through an averaging 
process similar to the one described above for the transverse momentum distribution of 
Drell-Yan pairs : one can factorize the b-distribution in eq. (|3~0"D , by evaluating Abn with 
q m ax at its mean value, i.e. write 

n(b, s) = n sof t(b, s) + A BN (b, < q max (s) >)o>t (31) 
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with 



and 



a 



jet 



zJ%f%*M*MW*£i (32) 



^ ( s ^_ y/i Sij J /i/a(xi) J - f) ,,,, 

W(SJ 2 J ^/, /a (xO J ^f j/b (x 2 ) f(dz) m 

with the lower limit of integration in the variable z given by z min = 4pf min / (sxix 2 ). Using 
for the valence quarks the same approximation as in the Drell-Yan case, we obtain to 
leading order 

3 fs 

< q ma x{s) >~ -PtmiJn- (34) 

for 2p t min « \/s. For p t min = 1-4 GeV, as in typical eikonal mini-jet models for proton- 
proton scattering |J, one obtains values of < q ma x( s ) > which range from 0.5 to 5 GeV 
for -y/i between 10 GeV and 14 TeV respectively A more precise evaluation of the above 
quantities depends upon the type of parton densities one uses, and will be discussed in a 
forthcoming paper. 

From the discussion about the large b-behaviour of the function h(b), we then expect 
ABN(b, s) to fall at large b more rapidly as the energy increases from yfs = 10 GeV into 
the TeV region. In Fig. 2, we compare this behaviour for the function A(b) with the one 
obtained through the Fourier transform of the squared e.m. form factor of the proton, i.e. 

^-I^^ 1 -^^ ^-O.llG^ (35) 

which is the one proposed by L.Durand et al. in the first eikonal mini-jet model for 
proton-proton collisions. The function A (b) from the Bloch-Nordsieck model is calculated 
for A = 0.1 GeV and values of < q max > which include those obtainable from eq. (j54]) in 
the energy range \/s ~ 10 GeV -r 14 TeV. 

We notice that, as the energy increases, A(b) from the form factor model remains 
substantially higher at large b than in the Bloch-Nordsieck case. As a result, for the same 
a j e t the Bloch-Nordsieck model will give smaller n(b, s) at large b than the form factor 
model and a softening effect of the total eikonal mini-jet cross-sections can be expected. 

In summary, a soft gluon summation model allows one to obtain a value for the intrinsic 
transverse momentum cut-off required in the b-distribution of partons in a proton in agree- 



ment with current phenomenological estimates [E2|. For such a description, an ansatz is 
necessary for the behaviour of the QCD coupling constant in the near zero region. By using 
suitable energy scale for the maximum energy allowed for emission of a single soft gluon 
by the valence quarks, and borrowing the expression for a s used elsewhere which leads to 
calculable non-singular integrals, we have obtained a progressively softer distribution in 
the large b-region as energy is increased. The physical interpretation of this effect can be 
traced to the fact that, as the energy increases, partons undertake scattering at smaller 
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Figure 2: Comparison between the b-distribution from the Bloch-Nordsieck model (full) 
and the form factor model (dots). 

and smaller b-distances. A comparison with the b-distribution from proton form-factor 
type models indicates a distinctly different behaviour in this large b-region suggesting a 
softening of the rise of the total cross-section in mini-jet models relative to the ones with 
the proton form-factor. Work is in progress to apply this model to hadron-hadron and 
photo-hadron collisions. 

This Work has been partially supported by CICYT, Contract # AEN 94-00936, EEC 
HCMP CT92-0026, and US-Department of Energy. 
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